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We study the coalescence of nanoscale metal clusters in an inert-gas atmosphere using const ant- 
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energy molecular dynamics. The coalescence proceeds via atomic diffusion with the release of 
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surface energy raising the temperature. If the temperature exceeds the melting point of the coa- 
lesced cluster, a molten droplet forms. If the temperature falls between the melting point of the 
larger cluster and those of the smaller clusters, a metastable molten droplet forms and freezes. 
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The ability to predict and control the size and structure of nanoparticles is of great 
importance for applications, as the arrangement of atoms is a key determinant of the optical, 
electronic and thermal properties of such particles. However, the structure of an atomic 
cluster often differs from that of the corresponding bulk material As the number of 
surface atoms is comparable to the number of interior atoms in a cluster, the surface energy 
plays an important role in determining the overall structure. The delicate balance between 
surface and internal energies often produces a complex dependence of equilibrium structure 
upon cluster size 

as. 

In an inert-gas aggregation (IGA) source, where metal clusters are produced by evapo- 
ration and then condensation from a vapor in an inert gas atmosphere, experiments have 
shown that the final cluster structure depends on the growth rate 4j. Molecular dynam- 
ics simulations of silver clusters support this finding, with studies of both the growth of 
clusters from small seeds ^| and the freezing of clusters from molten droplets finding a 
dependence of structure on kinetic factors. These examples illustrate that it is important 
to understand cluster kinetics if one wishes to predict cluster structure. 

The growth of clusters by coalescence is a kinetic process which will also effect structure. 
Indeed, the aggregation of atomic clusters has generated much interest due to the variety of 
different large-scale structures that can form, and the relevance of this to the functionality 
of devices built from clusters. For example, the deposition of clusters onto surfaces can 
result in highly ramified fractal structures p|, highly symmetric droplets j^, as well as 
many intermediate structures A recent proposal [H| uses cluster deposition to form 
wire-like cluster chains resembling nanowires with conduction of electrical current between 
neighboring particles. Understanding cluster coalescence is therefore crucial for reliable 
design and fabrication of such devices. 

Further, the distribution of cluster sizes produced in IGA sources suggests that coales- 

nn 

cence is an important growth mechanism pjj, I12j . Studies by electron diffraction of lead 



clusters grown in an IGA source have found structures that could not be fitted by standard 
model structures but which might be accounted for by coalesced structures [13]. To predict 
the structures of metal clusters produced in IGA sources we need to study the coalescence 
of free metal clusters. 

The collision and coalescence of small clusters, containing tens of atoms, has been studied 
extensively using molecular dynamics. For example, an early study showed that the decrease 
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in surface area during fusion can lead to an increase in kinetic energy 14| . Simulations of the 
coalescence of larger clusters have been predominantly conducted at constant temperature 
[15I [lfl, which probably best represents the coalescence of supported clusters with the 
substrate acting as a thermostat 16]. However, as we will argue below, the coalescence 
of nanometer-sized free clusters in an IGA source is best described by constant energy 
molecular dynamics. We are aware of only one such study, which considered the fusion of 
both pairs of liquid and pairs of " glassy" solid silicon clusters [l7 1 , and demonstrated that the 
kinetic energy generated during coalescence can considerably shorten coalescence times by 
enhancing atomic diffusion. Recent kinetic Monte Carlo simulations have shown, however, 
that the coalescence of facetted clusters can be considerably slower than surface-roughened 
particles jl^ . 

Our purpose is to study the coalescence of facetted solid lead clusters using molecular 
dynamics under conditions similar to those found in an IGA source. We examine the co- 
alescence of lead clusters using the empirical glue potential developed by Lim, Ong and 
Ercolessi Recently it was shown that many-body effects in the glue potential pro- 

duce unusual surface- reconstructed icosahedra in simulations of freezing [2Jj. Further, these 
surface-reconstructed icosahedra are thought to be globally optimum structures for the glue- 



potential for lead over a wide range of sizes [20, l2lL l22j . It is of interest then to study the 
coalescence of such icosahedra under conditions that might occur in an IGA source. 

An estimate of the mean time, r 9 , for a cluster of radius R between collisions with an 
inert-gas atmosphere of atomic mass m g , pressure P and temperature T in the source can 
be made as follows j^|: 

9 Pi? 2 V 87T • 1 } 

For a 2 nm radius lead cluster in a helium atmosphere at 5 mbar and 500 K, the time between 
collisions is approximately 1 ns. Thus over periods of a nanosecond or so, a constant-energy 
simulation of cluster coalescence will be a reasonable approximation to conditions in an IGA 
source. We will use constant-energy molecular dynamics here. 

Two 565-atom surface-reconstructed icosahedral clusters are initially prepared by con- 



struction as detailed in Hendy and Doye 



22j • Each of these clusters consist of a 147 Mackay 



icosahedral core with two reconstructed outer shells which differ from the usual Mackay 
icosahedral surface termination. Each cluster is then equilibrated in isolation during a 
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constant-temperature simulation for 10 5 time-steps at some common initial temperature Tj. 
The clusters are then placed in a common simulation cell so that the edges are aligned and 
the corresponding edge atoms are 0.5 nm apart (approximately next-nearest neighbor dis- 
tance but within the glue potential cut-off). In this work, we will not examine the effects of 
relative orientation, initial structure or impact parameter of the clusters in detail. Further 
simulations reveal that changing relative orientations does not alter the conclusions reached 
here, although we have not investigated the case of non-zero impact parameter. 

The constant-energy simulation then proceeds at a total energy equal to the sum of the 
final energies from each isolated equilibration run. Here we use a 3.75 fs time step which 
ensures energy is conserved to within 0.01 eV/atom over a 10 ns simulation. Figure Q shows 
a series of snapshots of the initial stages of coalescence between two 565-atom clusters each 
equilibrated at 300 K. Coalescence typically begins tens of picoseconds into the simulation 
when an atom on one of the aligned edges makes contact with the adjacent edge of the other 
cluster. The first snapshot in figure H shows the situation at this initial moment of contact, 
with the formation of an initial neck shown in subsequent frames. Neck growth is essentially 
complete by the final snapshot in figure ^ Figure El shows the evolution of temperature and 
the aspect ratio 2^ 

As the neck grows, the surface area of the merging clusters decreases, leading to the 
release of considerable surface energy. This energy is converted to heat, as can be seen in 
Figure El by the large increase in temperature during neck growth. After initial contact and 
the period of rapid neck growth, the aspect ratio evolves more slowly. In Figure we can see 
that the cluster is essentially spherical in shape after a period of 0.5 ns after initial contact, 
as the aspect ratio is very close to one by this stage. Examination of cluster structure reveals 
that it is static after this time. We conclude that coalescence is complete within 1 ns. The 
final structure of the 1130-atom cluster is that of a surface-reconstructed icosahedron (see 
Figure El). 

The amount of surface energy released as heat will clearly play a role in enhancing atomic 
diffusion during coalescence. In fact, if the initial temperature of the two clusters is close 
enough to the melting point, the final coalesced cluster may be molten. Such a case is shown 
in Figure El which illustrates the evolution of temperature and aspect ratio of two clusters 
initially at 420 K. Here the coalescence is far more rapid than in the previous example, 
occurring in tens of picoseconds rather than over 0.5 ns. The clusters form a neck and 
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FIG. 1: Coalescence of two 565-atom icosahedra initially at 300 K. The sequence of snapshots at 
3.75 ps intervals shows the early growth of the neck after the initial contact. 



then melt, with the coalesced cluster rapidly becoming spherical. Note the initial spike in 
temperature prior to melting, followed by a drop in temperature due to the heat of melting. 
The final cluster is molten at a temperature of approximately 480 K. In the IGA source, 
this coalesced cluster will cool slowly and freeze before reaching equilibrium with the inert 
gas atmosphere. 

It is instructive to consider simple models of coalescence and melting. Consider the 
coalescence of a spherical particle of radius R\ with a spherical particle of radius R2 > R\. 
Conservation of energy tells us that the increase in thermal energy of the coalescing cluster 
must balance the reduction in surface energy. Thus, the increase in temperature of the 
aggregate can be estimated as: 



where a is the surface tension, c v is the heat capacity and p is the density of the two particles. 
In general, a, c v and p will depend on cluster size. However, we will ignore this dependence 
here and use the bulk values for these parameters. For the simulation presented in figure |2l 
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FIG. 2: Evolution of the temperature and the aspect ratio during the coalescence of two 565-atom 
icosahedra at 300 K. At the point of first contact between clusters (approximately 30 ps into the 
simulation), the temperature rises rapidly as surface energy is released. The inset shows the final 
cluster structure. 

R\ — R'2 — -R(565) = 1.6 nm and equation (0) predicts AT = 160 K, which is close to the 
value attained in the simulation of 150 K. This increase in temperature is maximized when 
the ratio of change in surface area to volume ratio is maximized i.e. when Ri = 0.852i? 2 for 
a given i?2- If T + AT exceeds the cluster melting point then we would expect the cluster 
to melt during coalescence. 

However, it is well-known that the melting point, Tm(R), of a metal cluster depends 
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FIG. 3: Evolution of the temperature and aspect ratio during the coalescence of two 565-atom 
icosahedra initially at 420 K. At the point of first contact between clusters (approximately 30 ps 
into the simulation), the temperature rises rapidly as surface energy is released. The final state of 
the coalesced cluster is a molten droplet shown in the inset. 



strongly on the radius, R. A simple model for this dependence can be written as 



T M (R)=T B [l-^) (3) 



RJ 

where L is some length which depends on bulk material parameters. Obviously, this model 
fails for clusters of radius R ~ L and does not account for the irregular variation in melting 



point observed in small clusters [24J, but experimental results for lead [25], where L ~ 0.5 



nm, suggest it is adequate for R > 2 nm. By simulation we have identified the melting 



point of the initial 565-atom clusters to be Tm(565) = 400 ± 10 K and that of the coalesced 
1130-atom clusters to be Tm(1130) = 480 ± 10 K. Equation (jSJ) predicts the melting point 
for the 565-atom clusters (R = 1.6 nm) to be 410 K and for the 1130-atom clusters (R = 2.0 
nm) to be 450 K. This rough agreement between Q and simulation will be sufficient for our 
purposes here. 

Clearly one can always arrange an initial temperature T, such that T + AT > Tm{{R\ + 
R 3 ,) 1 ^ 3 ). However, for clusters of equal size R\ = R 2 , comparison of equation © and 
equation (jSJ) suggests that there is a critical cluster size R = 1 nm, below which coalescence 
will always cause melting (i.e. since at Ti = 0, AT > T M ((Rl + R 3 2 ) 1/3 ))- However, we have 
not tested this using simulations here as the glue potential (l9| may be unsuitable for such 
small clusters j^j]. 

Figure 0] shows the final temperatures of a series of pairs of 565-atom clusters that have 
undergone coalescence. The initial temperature of the pairs has been varied from 250 K to 
400 K. The solid line is the value of AT predicted by equation (j2j). We see from the figure, 
that equation (j2J is adequate (considering the approximations involved) until an initial 
temperature of T = 340i^. Above this initial temperature the final state of the coalesced 
clusters is molten. The highest final temperature reached in this series of simulations is 
Tj + AT = A78K. Thus when the melting point of the final cluster is exceeded during 
coalescence, a molten cluster results. The drop in the final temperature is due to the latent 
heat of melting. Note, that as the initial temperatures of the simulations correspond to the 
total energy, figure |U is essentially a caloric curve for the coalesced clusters. The coexistence 
of both liquid and solid clusters between final temperatures of 425 K and 475 K corresponds 
to the van der Waals loop in the microcanonical ensemble j^j]. 

In figure 03 we illustrate the coalescence of two 565-atom icosahedra initially prepared at 
330 K. Coalescence occurs rapidly, with the aspect ratio reaching one 0.2 ns after initial 
contact. There is a corresponding increase in temperature of approximately AT ~ 90K at 
this time. Inspection of the cluster structure at this time reveals that it is highly disordered 
and is consistent with that of a molten structure 2i3]. The smaller than expected increase in 
temperature during coalescence also suggests the cluster is molten. However this structure 
is evidently unstable as at t ~ 0.5 ns, the cluster structure undergoes a transition to form 
a more regular icosahedral structure. There is a corresponding increase in temperature at 
this time. This route to coalescence via an unstable molten intermediate state was found 
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FIG. 4: The final temperature of 75 pairs of coalescing 565-atom clusters as a function of initial 
temperature. The solid line shows the final temperatures predicted by equation (J2J). At initial 
temperatures of 340 K and above the final coalesced cluster is molten and the final temperatures 
are lower due to the latent heat of melting. 

to be typical of coalescing clusters at initial temperatures of 320-330 K. It is interesting 
to note that the temperature, T = 420 K, of the unstable droplet falls between Tm(565) 
and T M (1130). This suggests a crude interpretation of this process, namely that the two 
individual 565-atom clusters melt when T exceeds Tm(565) during coalescence, but that the 
resulting 1130 atom molten droplet is now unstable as T < T/u(1130). Finally, it seems 
reasonable to characterize the liquid state seen immediately after coalescence as metastable: 
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FIG. 5: Evolution of temperature and aspect ratio during the coalescence of two 565-atom icosahe- 
dra initially at 330 K. The initial coalesced structure is liquid and evidently unstable. The cluster 
subsequently undergoes a structural transition to form an icosahedral structure. The inset shows 
the final cluster structure. 

while for small clusters dynamical coexistence between the solid and liquid state has been 
observed J27 1, the energy barriers between the liquid and solid states scale as N 2 ^ 3 J^, 
hence repeated melting and freezing of the 1130-atom cluster will not be observed. 

To summarize, we have identified three regimes for the coalescence of free solid metal 
clusters: 1. solid-solid coalescence by diffusion. This will produce highly-defective clusters, 
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perhaps with large numbers of stacking faults (see Ref. [la]). 2. complete melting of solid 
clusters as they coalesce to form a single liquid droplet. In an IGA source, the final structures 
will depend on the cooling rate due to equilibration with the inert gas 3. complete 
melting of solid clusters coalescing to form an unstable liquid droplet which subsequently 
solidifies. The solid structures produced by this mechanism will resemble those produced by 
a rapid quench. 
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